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Abstract. The ac driven sine-Gordon equation is studied analytically and numerically, with the aim of 
providing a full description of how soliton solutions behave. To date, there is much controversy about 
when ac driven dc motion is possible. Our work shows that kink solitons exhibit dc or oscillatory motion 
depending on the relation between their initial velocity and the force parameters. Such motion is proven 
to be impossible in the presence of damping terms. For breathers, the force amplitude range for which 
they exist when dissipation is absent is found. All the analytical results are compared with numerical 
simulations, which in addition exhibit no dc motion at all for breathers, and an excellent agreement is 
found. In the conclusion, the generality of our results and connections to others systems for which a similar 
phenomenology may arise are discussed. 

PACS. 03.40.Kf Waves and wave propagation: general mathematical aspects - 74.50.+r Proximity effects, 
weak links, tunneling phenomena, and Josephson effects - 85.25.Cp Josephson devices 



1 Introduction 

In the past two decades, solitons have become a very 
generic and useful paradigm of intrinsically nonlinear phe- 
nomena in many different fields of physics Gener- 
ally speaking, solitons can be classified B0] in three main 
groups: dynamical, envelope, and topological solitons, ac- 
cording to their most relevant characteristics such as sta- 
bility, defining parameters, response to perturbations and 
so on. A canonical example of the most important of those 
groups H, topological and envelope solitons, in physical 
systems is the sine-Gordon (sG) equation, widely used 
both because it is integrable (meaning that the corre- 
sponding initial value problem can always be solved ana- 
lytically in closed form) and because it accurately models 
a variety of problems and applications. With respect to 
the integrability of the sG equation, this property gives 
us the possibility to know exactly its multisoliton solu- 
tions, that can be classified as kinks (topological) and 
breathers (envelope) . On the side of the applicability, the 
sG equation describes propagation of ultra-short optical 
pulses in resonant laser media Q , a unitary theory of ele- 
mentary particles [f7j,|sl, |9l,|l0|l , propagation of magnetic flux 
in Jose phs on junctions |11||, transmission of ferromagnetic 
waves p2| , epitaxial growth of thin films motion 
of dislocations in crystals |ll||l6||, flux- line unlocking in 
type II superconductors [|l7j, DNA dynamics [ jT^|ls| .pu||. 
and many others. 

When using the sG equation as a model for an actual 
physical situation, it is often necessary to account for fac- 
tors that cause deviation from the perfect system, arising, 
for example, from forces acting on it, thermal effects, fluc- 



tuations, dissipation, or spatial modulations (determinis- 
tic or random). To account for some or all of those, ap- 
propriately chosen perturbing terms have to be included 
in the sG equation to quote a few instances of such 

perturbations, (see also the next paragraph) let us men- 
tion the studies of forces acting over a DNA molecule, or 
long DNA fragments containing regions of finite size and 
specific structure [^l|, additive and multiplicative noise 
sources fl22| , spatially periodic parametric potential [^3| 
or damping with spatiotemporal periodic driving J24| . The 
work reported on here belongs to this class of problems; 
specifically, it focuses on the action of ac (sinusoidal in 
time, homogeneous in space) forces on the solitonic solu- 
tions, kinks and breathers, of the sG equation. In mathe- 
matical terms, this means we have to deal with the equa- 
tion 

<ptt-<l>x X +sm(4>) = -l3<j> t +f(t), (1) 

where subindices indicate derivatives with respect to the 
corresponding variables, P4>t is the usual damping term 
and f(t) = esin((5i + Sq) is an external periodic force de- 
scribing, for example, a long Josephson junction under the 
application of a uniform microwave field [lljj . From now 
on, the amplitude of the force, e, will be our perturba- 
tion parameter, i.e., it is assumed that the force is not too 
large. We will see, however, that we can deal with reason- 
able and useful values of e, so our calculations and results 
are not purely academic or mathematical but rather they 
have physical relevance. 

The studies of forced sG systems date back to the late 
seventies, to the seminal works of Fogel et al. J2j| and of 
McLaughlin and Scott Eq] where the ideas of collective 
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coordinate techniques (see |2| for a review on those) were 
first introduced to deal with force terms added to the sG 
equation. Their results, confirmed by numerical simula- 
tions of the full partial differential equation, showed that 
dc (constant both in space and time) forces act on the sG 
solitons as if they were point-like particles, accelerating 
them up to the maximum velocity allowed by the equation 
in the absence of damping or to a terminal velocity in the 
presence of damping. The problem of dc forces was given 
further consideration in subsequent works, such as p7|,pc|, 
where aspects such as breather instabilities and pattern 
formation were the main subjects; as our principal goal 
is the study of ac forces, we will not go into more detail 
on those here. Dissipation effects on breathers were first 
studied in |29[] , where McLaughlin and Overman showed 
that free (i.e, not driven) damped breathers where anihi- 
lated, their energy being dissipated into radiation. Sub- 
sequently, Lomdahl and Samuelsen |3^,|l| showed that 
the introduction of ac forces in the damped sG equation 
compensated those losses and stabilized breathers, whose 
frequency became modulated by the driving. This is the 
main result about ac breather dynamics, and if fact we 
are not aware of later relevant work on ac driven breather 
dynamics. On the other hand, the kink case has a longer 
and more problematic story, that begins with the research 
by Olsen and Samuelsen in |3^], where the authors stud- 
ied the effect of dc and ac forces over sG kinks. In their 
analytical calculation, they considered two special cases 
for an ac force given by esm(5t + 5q), namely 5q — 0, tt/2 
with (5^1 and sufficiently small kink velocities. Their 
results were obtained in the approach that the solution 
of the perturbed sG equation can be divided into a kink 
part and a vacuum part. Although most of the paper deals 
with dc forces, the authors conclude that for the two cases 
studied for ac forces, one leads to kink dc motion (So = 0) 
whereas for the other (So = tt/2) only oscillatory motion 
is found. For both choices, the work is further restricted 
to null initial velocity only, u(0) — 0, with the subsequent 
loss of generality of their findings. Numerical simulations 
confirmed the predictions and the accuracy of the analyti- 
cal calculations. However, no other initial phase differences 
were studied, and therefore a general analysis (including 
also initial velocity) was lacking. In addition, they did not 
present any results for ac driving plus damping, saying 
simply that similar results could be obtained in that case. 
A few years after that work, Bonilla and Malomed |£| 
revisited the same problem and concluded that dc motion 
of ac driving kinks was only possible in discrete sG equa- 
tions (i.e., the spatial variable is discretized), and only in 
the presence of damping. Such a result, in contradiction 
with that in (32) , was later shown to originate in an incor- 
rect analytical approximation to the problem in [ p4| ; the 
perturbative theory worked out in that paper suggested 
that kink dc motion induced by ac forces in the presence 
of damping was not possible, and the authors stated their 
belief hat such phenomenon was unlikely in general, i.e., 
in the undamped case (at that time, they were not aware 
of the results of Olsen and Samuelsen p2[). Finally, the 
last contribution to this topic prior to our work appeared 



in [ p5[ , where dc motion of undamped ac driven kinks was 
found in simulations of the Frenkel-Kontorova (discrete 
sine-Gordon) model. In this case, the phenomenon was 
explained in terms of a resonance between the external 
force and internal oscillatory modes of the discrete kink 
(absent in the continuum models). 



The brief summary in the above paragraph is enough 
to understand that the knowledge on the ac driven dynam- 
ics was not complete and systematic and, what is worse, it 
severely suffered from the contradictory results obtained 
about the possibility of dc propagation of kinks. In view 
of this, we first addressed this issue in J36|, a prelimi- 
nary work on small (non-relativistic, \u(t)\ <C 1) velocity, 
ac driven kinks. Using the techniques developed in p6| , 
we considered the general ac force esin(<5i + <5o) with and 
without dissipation. By this means, we were able to show 
analytically that there indeed exists dc motion of kinks 
provided dissipation is absent. Moreover, we showed that 
the velocity of such dc motion depends on the relation 
between initial velocity and ac force parameters (ampli- 
tude, frequency and phase), thus generalizing the results 
of Olsen and Samuelsen 32 for the non dissipative case. In 
the presence of damping, we also showed within the non- 
relativistic approximation that kink dc motion is never 
possible. This conclusion is different from that of Olsen 
and Samuelsen |3^| , who said "similar [to the undamped 
case] results are obtained", not a very explicit statement; 
it is also contrary to the claims of Bonilla and Malomed 
p3[ , and finally it confirms the results of Cai et al. p4j . 
Numerical simulations fully reproduced the predictions of 
the approximate calculations establishing their validity to 
a high degree of accuracy. In the present paper we consid- 
erably extend the study reported in |36 in a systematic 
way, by analyzing i) the relativistic case for kinks, and ii) 
the effects of pure ac driving on breathers, that had not 
been studied so far. Our main conclusion will be that the 
application of external ac driving forces can cause the dc 
motion of sG kinks (for the whole range of initial veloci- 
ties) depending on their phase; for breathers, we show that 
in the absence of damping they can be self-sustained up 
to a certain threshold of the force intensity, at which they 
break up into kink-antikink pairs. We also prove that com- 
bining the damping effect with the ac force the dc motion 
of sG kinks becomes asymptotically an oscillatory one. In 
Sec. |[ we use a collective coordinate approach to obtain 
the equations of motion for the collective variables finding 
that the mean velocity of the kink (antikink) and the fre- 
quency of the breather are functions of the parameters of 
the ac force. We compare our analytical results by means 
of numerical simulations in Sec. S and verify that they 
are in an excellent agreement. Finally, Sec. |i| summarizes 
our findings and conclusions, and presents a brief discus- 
sion about some other systems in which the phenomenon 
described in this paper can also appear. The paper closes 
with two appendices, one devoted to the proof of a mathe- 
matical intermediate result needed in Sec. and the other 
analyzing the effect of another dissipative term (with con- 
stant driving), u xxt , of particular significance for Joseph- 
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son junction dynamics, to show the power of our technique 
(see Sec. ^ to solve the relativistic problem. 



2 Collective coordinate approach 
2.1 Kink dynamics 

We begin our study by addressing the problem of kink 
dynamics, considering Eq. (]l|) when there is a single kink 
present in the system with initial velocity u(0), t = be- 
ing taken as the time at which the ac force is switched on. 
In order to study the subsequent kink evolution governed 
by Eq. ([!]), we apply the collective coordinate perturba- 
tion theory in the version introduced by McLaughlin and 
Scott |26f| . If j3 and e are small parameters we may assume 
that the solution of Eq. (|l|) has the same form as the un- 
perturbed sG equation [Eq. (Q) with (3 — e = 0], except 
that now we allow x (t),Xo(t) and u(t) to be functions of 
time, so that 



(j)(x, t) = 4 arctan exp ± 



X (t) = / u(t')dt' 
Jo 



X - X (t) - Xg(t) 




where the positive (negative) sign corresponds to a kink 
(anti-kink) solution. We note that X(t) = x (t) + Xo(t) 
[u(t)] has the meaning of the position (velocity) of the cen- 
ter of the soliton, and that the main assumption under- 
lying this approximation is that radiation effects induced 
by the perturbation are neglected. This will be verified a 
posteriori by comparing with the numerical simulations. 

With this Ansatz, we proceed much in the same way 
as McLaughlin and Scott |^6|: we first compute the vari- 
ation of the energy and momentum of the unperturbed 
sine-Gordon system due to the perturbation, and applying 
them to Eqs. ^ and (|J) obtain those variations in terms 
of our unknown functions. Imposing that it is consistent 
with the result of computing the variation of energy and 
momentum from the kink expression with constant pa- 
rameters and making them functions after the calculation 
(see [|6| for details; see also Q for other versions of the 
collective coordinate technique) , we arrive at the following 
equations for the position and velocity of the perturbed 
kinks (the ± sign corresponds to a kink or an anti-kink, 
respectively) : 

^ = -1(1- u 2 )tWl - u*ef(t) + 4/3t»], (4) 



u(t = 0) = u(0); 
dX 

If = U ^ 
X(t = 0) = X(0). 



(5) 



(equivalently, Eq. ([|) could written as x (t) — const). For 
the dissipation-free case [f3 = in Eqs. @) and (Eh], Eq. 



(||) can be exactly solved, yielding 

u{t) = F[u(0),5 ]/ (1 + F[u(0), S } 2 ) 1/2 , (6) 
«(0) 



F[u(0),6o] 



Vi-« 2 (o) 



± 



± — [cos(<Ji + Jo)-cos(i )]- (?) 



Once we have obtained the velocity, we need to solve Eq. 
(||) to find the kink trajectory. Unfortunately, it is not pos- 
sible to find an analytical solution, so in order to get some 
insight about it we restrict ourselves to the non-relativistic 
limit, i. e., \u(0)\ <C 1, and imposing |7re/4<5| <C 1, (which 
is no loss of generality as we have already assumed f«l) 
we can find an approximate solution for the position of 
the center of the kink, X(t), that reads 



X{t) ~ X(0) 

7re[sin(<5i 



7TC 



u(p) T 77 cos(<J )] t ± 
4o 



±- 



S ) - sm(5 )} 



AS 2 



(8) 



Equation (Bh tells us what the kink trajectory is, at 
least while its velocity is not large. Keeping in mind this 
caveat, it is evident from the above expression that only 
for the "resonant" velocity 



u(0) 



±[^cos(<5 )], 



(9) 



the linear term will vanish and the motion of the kink (or 
antikink) will be oscillatory; for any other velocity values 
(fixed 5, 5q and e) the ac driving force induces a dc mo- 
tion of sG kinks, even for the kinks initially at rest. This 
result contains those found by Olsen and Samuelsen [ |32| 
(also obtained in the non-relativistic approximation, by 
the way) as particular cases for two choices of 5q. 

So far, we have developed a collective coordinate the- 
ory in the non-relativistic limit. This have already pro- 
vided information on the kink dynamics, but, in addition, 
we can benefit from that to find the unrestricted relativis- 
tic result, by virtue of the following fact: It can be proven 
(see Appendix I) from equations M) and (0) that the re- 
lation 



F[u(0),6 ]dt = 0, 



(10) 



is a necessary and sufficient condition for the motion to 
be oscillatory with period T = 2w/S. Hence, integrating 
( |Io|) we find that kinks will oscillate if and only if 



7= 5?=; = ± 77 c os(<5o)]- 
v /l-u 2 (0) 4d 



(11) 



Let us now discuss this condition. To begin with, the 
condition ( pd| ) involves four parameters, so we can fix, say, 
5, e and u(0) and then use ( fj"l| ) to find the value of 5q 
such that the motion becomes oscillatory. Second, the only 
change with respect to the approximate condition (^|) is 
the appearance of the Lorentz factor 7 — yj\ — u(0) 2 in 
Eq. (|l l|) , which vanishes in the small velocity limit, i. e., 
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the result (|^) is recovered. Equation is very impor- 
tant, because it implies that (within the collective coordi- 
nate approach) , any ac force applied on a kink will induce 
its dc motion with a velocity directly related to their rel- 
ative phase [vanishing only for a very special choice of 
that parameter, So]- Our perturbative calculation coin- 
cides with that of Olsen and Samuelsen in the limit and 
special cases studied by them, thus generalizing their re- 
sults in a comprehensive manner. We postpone the nu- 
merical validation of our findings to Sec. |3|. 

After dealing with the undamped case and arriving 
to an expression valid for all kink velocities, let us now 
turn to the damped problem, including a damping term 
(P 7^ 0) in the ac driven sG equation. We will now demon- 
strate that dc motion of sG kinks is not possible in this 
case. Indeed, in our previous work p6[ we analizcd the 
effects of damping {(3 ^ 0) and ac force non-relativistic 
kinks and found that kinks exhibit dc motion except for a 
transient, after which they reach a final oscillatory state 
around a point whose location depends on the initial con- 
ditions. However, nothing was said there about the general 
case, with no limits on the initial velocity, so we need to 
see what is the behavior of relativistic kinks by looking 
back again at Eq. (|J). The result found above for the un- 
damped case suggests to use y(t) = u(t)/y/l — u(t) 2 as a 
new variable; by doing so, Eq. (|J) can be linearized and 
the analytical solution for the velocity computed. We thus 
arrive at 



u(t) 
r(t) 



•(f) 



c exp(—(3t) =F 



(12) 



4(/3 2 + S 2 ) 

x[/3sm(6t + 6o) - 6 coa(5t + 5 )], (13) 

ti(0) £7T 
C = y/l - u(Q) 2 W+P) X 

x[/?sia(£o)-£cos(£ )]. (14) 

From Eqs. it can be easily seen that the asymp- 

totic behavior of u(t) is given by 

m 



f(t) = ± 



ire 



(15) 



4(/3 2 + <5 2 ) 
x [8 cos(<5i + So) 



(3sm(8t + 8 )}. (16) 

Even if we cannot explicitly solve Eq. (|5|) when u(t) is 
given by Eqs. (12|, ^3p4| ) for the trajectory X(t), the above 
findings are all we need to prove (always within the col- 
lective coordinate formalism) our previous claim that kink 
dc motion is asymptotically forbidden. To this end, it is 
enough to show that for to and t large enough, X(t) ~ 
X(t) = X(to) + / t * u(t')dt', is a periodic function. This is 
equivalent to proving that 



t+T 



u(t')dt' = 0; 



A straightforward calculation shows that the condition 
( |l7| ) is always true, giving grounds to our conclusion that 
kinks can not exhibit dc motion in the damped, ac driven 
sG equation. This is yet another proof that the arguments 
of Bonilla and Malomed (33| supporting the opposite con- 
jecture were not correct. 



2.2 Breather dynamics 

Following our study for kinks, we have carried out an anal- 
ysis of the effect of ac driving over the breather solutions, 
based again on the collective coordinate approach (see de- 
tails in Ull)- 

We concern ourselves mostly with the 
undamped case, as the damped problem has been already 
carefully discussed in [|o[|3^]. The collective coordinate 
Ansatz wc will use for breathers is 



4>{x, t) = 4 arctan 



g(t) 



cosh[a; to(t)] 



(18) 



where the breather frequency is J?(i) = \J\ — lu 2 (t) . This 
Ansatz is assumed to be valid when (3 and e are small. 
In this case the collective variables are the two functions 
u>(t) and g(t), which can be found solving a system of 
differential equations equivalent to that in Eqs. ([|, ||) for 
the kinks. The procedure to find the equations is the same 
and yields 



~dt 



fit) 



4 v/l + .9 2 W 



x cos |ar(Mii 



w(t) 



9(t)} 



(19) 



U){t = 0) = LUq 

dg 
dt 



w{t) 



4uj(t) 



/(t)arcsinh[g(i)] \ x 



. ,Jl - UJ 2 {t) . s .. 
x cos(arcsm[- — g{t)\), 



(20) 



j(t - 0) = 0. 



We have not been able to solve the above equations. Never- 
theless, they are still useful in order to predict the breather 
stability conditions, as we now show. If u>(t) ^ 0, the sys- 
tem of Eqs. (H|,|^) leads us to 



/(*) 



1 diO 7T£ 

u~dt ~ 4^J 2 y/l +g*(t) 
1 



dg 



(1 



4w 2 



/(t)arcsinh[ 3 (i)]) dt 



(21) 



Straightforward algebraic manipulations of this equation 
yield 

f t H"(t)\ = |ln[l + i ^/Warcsinh( 5 [t])]. (22) 



where terms of order e 2 and eS have been neglected. The 
\ > advantage of this expression is that we can now integrate 
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it; integration and subsequent expansion of u)(t) around 
wq in powers of e allow us to write (up to orders e 2 again) 



0(x,t) 



w(t) =wo[l + <r(t)], 



7T6 

cr(t) = — 2 /(i)arcsinh 



u sin(- v /l -Ldpt) 



(23) 
(24) 



In order to proceed, we have to recall the following 
property: breather solutions can only exist provided that 
< u>(t) < 1 (which implies that the breather energy is 
bounded above and below, < Ht, re (t) < 16, kinks hav- 
ing energy equal to or larger than 8 in our dimensionless 
units). Having this in mind, we can find a condition for 
the transformation of a breather into a kk pair, that turns 
out to be 



Sthr 
tthr 
Cthr 



min (ei, £2) , if M > and m < 0, 
e 2 , ifM<0, 
ei if m > 0, 



with 



<U 
£2 



4cj (1 - u} ) 
7rm 



(25) 

(26) 
(27) 



m and M being the minimum and maximum respectively 
of the function a(t) [see Eq. above]. 

The physical meaning of the somewhat cumbersome 
result is the following: when sG breathers are perturbed by 
an ac force there always exists a value of e, e t ^, such that 
perturbations smaller than the threshold allow breathers 
to exist with modulated frequency fl(t) = yl — ^ 2 {t), 
where uj(t) is the function (|23|); on the other hand, when 
e > tthr the breather transforms into a kink-antikink pair 
or decomposes into a radiation, depending of the param- 
eters of ac force and its frequency. Therefore, we arrive 
at our major conclusion regarding (undamped) breather 
dynamics, very different from the kink analysis in the pre- 
vious subsection: Breathers are only stable for a limited 
range of force intensities, and in that range the effect of 
the driving is to modulate their frequency as discussed 
above. As before, we will have to compare this prediction 
to numerical simulations to verify its validity. 



3 Numerical results 

Although the collective coordinate approach has been al- 
ready tested very many times Q and is by now a standard 
technique to deal with soliton-bearing equations, the an- 
alytical expressions and the conclusions we have drawn 
from them should not be taken as valid before compari- 
son with detailed and accurate numerical simulation of the 
full partial differential equation (jl|). To this end, we have 
numerically integrated it by means of a standard fourth 
order Runge-Kutta algorithm |3q] with initial conditions 



*(x.t) 



*(x,t) 




Fig. 1. Verification of collective coordinate predictions in the 
absence of dissipation. Simulation starts from a kink with ini- 
tial velocity m(0) = 0.5 located at X(0) = 0, and subject to an 
ac force given by 0.1sin(0.1i + So), (a) 80 = 0.75, (b) So = 0; 
notice that the direction of motion is opposite in both cases, 
(c) So = 0.6, critical value exhibiting oscillatory motion. 



given by an unperturbed sine-Gordon kink or breather 
with unrestricted initial velocity m(0), and boundary con- 
ditions (f> x (L — ±50, t) = 0. We have carefully checked 
that our results did not depend on the choice of spatial 
and temporal steps, or the simulated system size (which 
can affect the simulations through the generation of spure- 
ous radiation). 

We first address the validation of the results for kinks, 
beginning with the undamped case. In Fig. [I] we show the 
time evolution of a sG kink with initial velocity u(0) — 
0.5 (which is certainly far from being small) for different 
phases of the ac force 5p and e = 0.1, 5 = 0.1. From Eq. 
( |IT1) we know that the kink (with initial velocity u(0) = 
0.5) will oscillate only if 5q — 0.745. Numerically, we 
searched for the critical phase of driving force and found 
Sq — 0.6 (see Fig. |l|c), i.e., the accuracy between these 
results it is of order of 0.1. This is a very good predic- 
tion if we consider the value of the velocity and the value 
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Fig. 2. Time evolution of the velocity of kink center mass. 
Parameter values as in Fig. 1. (a), (b) and (c) respectively. 
The solid curve shows the value of velocity, predicted by the 
collective coordinate perturbation [Eqs. (MM)], and the dotted 
line is the numerical integration of [Eq. (Pp] 



of e (indeed, the accuracy is of the order of e as should 
be expected). Note that there is practically no radiation 
visible in the simulations, which is in agreement with the 
correctness of the use of collective coordinate techniques. 
In connection with this, another interesting remark is that 
background motion, i.e., motion of the wings of the soliton 
induced by the ac force, is not visible in the simulations; 
this makes our approach very suitable for this problem, 
and as accurate as that of Olsen and Samuelscn (3^] who 
introduced a term accounting for this negligible effect. We 
note that smaller initial velocities compare even better 
with our theoretical results. To appreciate to a larger ex- 
tent the accuracy of our approach, Fig. || compares the 
analytical prediction for the velocity values [Eqs. (§,0)] 
with the velocity of the kink center of mass, obtained by 
integrating numerically Eq. ([l]). In these figures the pa- 
rameters of the perturbed sG equation were chosen as in 
Fig. |l[ We see from Fig. ||that the agreement is excellent, 
in spite of the fact that the prediction of the critical veloc- 
ity is not correct by a factor e. This is due, on one hand, 
to the fact that slight deviations from simmetry around 
u(t) — have dramatic consequences in the kink motion 
(note, e. g., Fig. ||c, the plot for the critical velocity, where 
the small discrepancy between the theory and the simu- 
lation is always in the same direction, towards negative 
u, corresponding to the theoretically predicted dc motion 
to the left), and, on the other hand, to the difficulties in- 
volved in numerically determining the critical velocity by 
checking that the kink just oscillates; this, among other 
things, depends on the integration time and on the time 
step. Hence, we believe that the correctness of our com- 
putations is much better appreciated by looking at the 
function u(t) than simply verifying the critical velocity 
value. 

We cannot close our study about sG kinks without 
providing numerical evidence that, as predicted by our 
analytical calculations above, it is not possible to induce 
dc motion on sG kinks by means of an ac force in the 
presence of damping. Figure |3| shows an example: a kink, 
initially moving with velocity u(0) = 0.5, perturbed by ac 
force 0.01sin(0.1t + 0.1) and a damping term (3 — 0.05, 
is stopped and ends up describing small amplitude oscil- 
lations around a finite region in space. This is but one 
example, although we must say that we have verified that 
our prediction is correct in very many other cases, with 
different choices of parameters, finding always an excellent 
agreement between analytics and numerics. As in the un- 
damped case, Fig. ||, depicts the predicted and simulated 
motions of the kink center, exhibiting a perfect agreement 
between both. It is of course reasonable that the collective 
coordinate approach works better for the damped case, 
because if any radiation arises in the system, the damp- 
ing wipes it away, leaving the kink as the only excitation 
present in the system, as the approximation requires. 

Finally, for the breather case we present in Fig. || a 
comparison between the analytical and numerical thresh- 
old values of e as a function of the initial frequency of the 
breather J? - We consider only the undamped case, (3 = 0. 
We have observed in our simulations the phenomenon pre- 
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0(x,i 




Fig. 3. Verification of collective coordinate predictions with 
dissipation. This is one example of soliton motion, initially 
with velocity tt(0) = 0.5, and subject to an ac force given by 
0.01 sin(0.1t + S ), with S = 0.1 and = 0.05. 

dieted by our analytical calculations: If the amplitude of 
the ac force exceeds the threshold value of e, the initial 
breather decomposes into a kk pair. As depicted in Fig. 

the collective coordinate result predicts very well the 
threshold amplitude in all the range of breather frequen- 
cies, even for rather large values of such amplitude, where 
the perturbative approach might not work in principle. 
We note that, for driving frequencies close to 1 (the lower 
edge of the phonon band of the sG system), the breather 
can as well decay into radiation, but we have not pur- 
sued this process in detail in view of the difficulties in- 
volved in distinguishing breathers from linear radiation in 
this limit. Therefore, we concern ourselves with the decay 
into kk pairs, this being the reason we have not gone over 
Qq = 0.9 [[39] in our numerical tests. Figure || shows an 
instance of our numerical determination of the threshold, 
for a breather with Qq = 0.2. We can see that if e = 0.0148 
(Fig. |^a) the frequency of the breather is modulated, but 
it remains a breather, as the oscillations of the field value 
at x = clearly show (of course, we have looked at the 
simulation over the whole spatial interval to confirm this). 
On the contrary, if we increase e (for example, e = 0.0149) 
then the breather breaks up into a kk pair (Fig. |6|b), a 
phenomenon which reveals itself in the suppression of the 
oscillations at x = indicating that the breather stops its 
periodic motion. We have thus verified the validity of the 
collective coordinate results for the breather case as well. 



0.8 
1 0-6 



CD 




-0.2 1 1 

100 200 300 400 
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Fig. 4. Time evolution of the velocity of kink center mass. 
Parameter values as in Fig. 3. The solid curve is the analyt- 
ical prediction [Eqs. jl^ , |l^ , p^] and the dotted curve is the 
numerical integration of [EqlTlh]. 
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0.6 




Fig. 5. Illustration of the threshold value of e vs the initial 
frequency J?o. Solid line represent the values given by the per- 
turbation analysis and the points are the values, obtained by 
numerical simulation of perturbed sine-Gordon equation. The 
frequency and the phase of ac force are 5 = 0.1, So = respec- 
tively. 



4 Conclusions 

In this paper, we have shown both by analytical (collec- 
tive coordinate theory) and numerical means the following 
three facts: i) Kink dc motion is possible in the sG system 
when forced by a pure ac force, except for a very special 
choice of the initial velocity related to the driving phase 
by equation ([ll]); in this case, the kink center oscillates 
around its initial position; ii) kink dc motion induced by 
pure ac driving is never possible in the presence of damp- 
ing, and iii) the undamped sG system forced by an ac 
driving allows the existence of breathers, whose frequency 
becomes modulated by the external one, up to a certain 



critical value of the amplitude that depends on the driv- 
ing frequency and on the breather parameters. Above that 
amplitude, the breather decays to a kink-antikink pair. 

As already stated, this findings arise from a perturba- 
tive calculation of the collective coordinate type, whose 
accuracy for very many soliton problems is well estab- 
lished [pj. In this respect, we want to emphasize our main 
achievement in this study: We have proposed a change 
of variables which allows to linearize and fully solve the 
differential equation for the kink velocity without any re- 
striction. We note that in previous works |32|,[56| only the 
non-relativistic limit could be dealt with in an approxi- 
mate manner. It is clear that this analytical procedure can 
be useful in similar equations, an example of this being our 
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Fig. 6. (a) The breather with modulated frequency for the 
value of e = 0.0148.(b) If e = 0.0149 the breather break into 
kk pair. The others parameter values are i?o = 0.2, 8 = 0.1, 
5q — 0. Compare with the analytical threshold e (e t h = 0.0135), 
obtained in Eqs. 



treatment of the long Josephson junction equation, with 
a more realistic dissipative term and a constant driving, 
in Appendix II. Therefore, having in mind the accuracy 
of the collective coordinate approach, the fact that there 
were no further approximations involved in our calcula- 
tions, and the numerical confirmation of our predictions, 
we believe that our study settles down once and for all 
the question of kink motion in ac driven soliton bearing 
systems, particularly so for the sG equation. 

From these results, which clarify the phenomenology of 
the ac driven sG equation, we can try to draw some more 
general conclusions on soliton-bearing problems perturbed 
by an ac force. The findings on discrete sG systems |35| 
suggest that for our conclusions to apply in other, different 
systems, it is necessary that the soliton or solitary wave 
under consideration does not have one or several inter- 
nal oscillation modes which can be excited by the driving 
force. To check for the possible new phenomena or mod- 
ifications of our conclusions coming from the existence of 
those modes, we have carried out some preliminary re- 



search on the ac driven </> model. Generally speaking, the 
same conclusions hold, i. e., in the absence of dissipation 
(j3 — 0) the kink (antikink) will exhibit oscillatory motion 

U(0) — ' 36 : cos (So)], (28) 



Vi-« 2 (o) 



= ± 



V2<5 



whereas when dissipation is present we again find that 
kink dc motion is not possible. However, we have found 
numerical evidence that for values of the driving frequency 
close to the internal mode one, the kink behavior is much 
more complicated, even chaotic. We are currently working 
on understanding this effect analytically, and tentatively 
we relate it to the resonance of the driving and the in- 
ternal mode f40| . Similar resonance phenomena were ob- 
served for the sG kink in a harmonic potential well, when 
the driving had a frequency close to the natural one of the 
kink in the well [El], which reinforces our interpretation 
and expectations. This leads us to the conjecture that the 
results we have obtained for the sG problem can be true 
in general, provided as already said that there are no res- 
onances with other modes intrinsic to the excitation. Of 
course, there also remains the question of drivings with 
frequencies above the phonon lower edge. We believe that 
our conclusions will still hold in that situation, in so far 
the presence of the external force does not cause an inter- 
action between solitonic and radiation modes; it is to be 
expected that our calculations will be very accurate in the 
presence of dissipation, which will damp away any excited 
radiation. Finally, as shown in Appendix II below, we want 
to stress that the predictions and mathematical results we 
have obtained in this work have direct application in many 
physical systems, and among them in Josephson devices 
[more realistically described by Eq. ( p^) below]. It would 
be most interesting to verify the possibility of dc motion 
of solitons induced by pure ac driving in real systems like 
that, more so in view of the potential applications such 
a rectification effect might have. Another implication of 
our results has to do with soliton generation in systems 
described by sG equations: Indeed, by driving the system 
with an ac force strong enough, we could break up ther- 
mally created breathers into kink-antikink pairs, which 
under the influence of the driving would separate moving 
towards opposite directions. This would be a very simple 
and clean way of creating solitons in, e.g., Josephson de- 
vices. It thus becomes clear that experimental work pursu- 
ing this and related questions, which is certainly amenable 
with the present capabilities, would be most useful in as- 
sessing the importance of our work. 

We thank Renato Alvarez- Nodarse, Jose Cuesta, Esteban Moro 
and Franz Mertens for conversations on this work. Work at 
CISC (Leganes) has been supported by CICyT (Spain) grant 
MAT95-0325 and DGES (Spain) grant PB96-0119. 



Appendix I 

In this appendix we prove the necessary and sufficient con- 
dition for the sG soliton motion to be purely oscillatory in 
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the presence of undamped ac driving. In the body of the 
paper we stated that this is so if and only if the condition 
(|l0|), or equivalently, (|TT|) , takes place. 

On the one hand, Substituting ( fjl| ) in ([?]) and after- 
wards in (^) we obtain for the velocity function u(t) the 
following expression: 



u(t) 



±ff cos(6t + go) 
1+ [ff cos(5t + 5 )]' 



(29) 



which can be integrated, yielding the position of the kink 
center of mass X{t): 



1 fire sin((5t + Jn) 



(30) 



Hence, it is evident from the above equation that X(t) is 
a periodic function with period T = 2ir/8. 

On the contrary, suppose that the kink center of mass 



X(t) = / u(t')dt' 
Jo 



(31) 



oscillates with period T = 2tt/<5, i. e., X(t) = X(t + T); 
this can be put as 



1(a) 



where 



t+T 



a + bcos(St' + S ) 
t v /l + (a + bcos(St' + S ))' 



-At' = 0, (32) 



u(0) 7recos((5o) 



^l-u(0f 



id 



(33) 



and b = tte/46. Noticing that, first 1(0) = 0; second, if 
a < -\b\ 1(a) < and if a > \b\ 1(a) > 0, and third, 
I' (a) > for all values of a, we arrive at the desired result 
that 1(a) = if and only if the relation (pd|) [a = 0] holds. 



Appendix II 

In |^6| realistic model of a long Josephson junction is pro- 
posed, given by the following perturbed sG equation 



4>tt - 4>xx + sin(</>) = -04> t + a(f> xxt - 7, 



(34) 



where the term (34> t represents the dissipation due tun- 
neling of normal electrons across the barrier, a(j) xxt is the 
dissipation caused by flow of normal electrons parallel to 
the barrier and 7 is a distributed bias current. In this ap- 
pendix we apply our theory to this problem in order to 
show that our technique to solve the equation for the col- 
lective velocity can be applied to other problems. In addi- 
tion, as the equation above is realized very approximately 
by actual Josephson devices, by studying it we will pro- 
viding experimental means of checking our results as well 
as predicting behavior which might have an application of 
interest in that context. 



To begin with, we note that if a = this equation 
coincide with (|l|) letting 5 — 0. Therefore, substituting 
6 = and 7 = -esin(6 ) in Eqs. (|T|,[iJ]JlJ) we obtain 



u(t) 



do exp(— (3t) + d 
■s/1 + (d exp(-/3i) + df 



(35) 



where d = u(0)/y/l - u(0) 2 -d and d = ±7r 7 /4/3. When 
t goes to infinity, we have u(oo) — > d/Vl + d 2 . Notice that 
this expression coincides with the equilibrium solution Uoo , 
obtained in |2(|. 

For a^0 the same procedure as above leads us to a 
first-order ordinary differential equation for u(t), which is 
given by 

f t = ±n ( i _ ^3/2 _ M1 _ u 2 } i au _ (36) 

If we now introduce the same change of variable as in Sec. 
[| y (t) = u(t)/y/l - u(t) 2 , Eq. (11) becomes 



dt 



± 



IT"/ 



Integrating 



p + xf 



4 

) we find 
ln( Cl 



C9 + f)w 



y-xi 



(37) 



2 , 



where 



3xi 
2k[ 



3p 
2q 

D : 

Xi ■ 



arctan , 



a 

37T7 

=F - 

4a 

(VD-q) 1 ^ 3 



2 




- C 2 



= —at, 



(VD-qy/f 



k\ = 3p ■ 



Cl 



(1/(0) 



^) 2 + k 2 



C2 = arctan 



2/(0) - x x 
2/(0) + 



and 



fei 



(38) 

(39) 

(40) 
(41) 
(42) 

(43) 

(44) 
(45) 



For large values of t, 2/(00) — > x\ , so the kink will move 



with constant velocity u(oo) 



Zl , which depends on 



the parameters a, f3 and 7. For example, if 7 = (q = 0), 
w(i)|t->oo = 0, i.e. by the influence of both dissipation 
terms (—f34>t and a(j> xxt ) the kink will be stopped. 



References 

1. M. Remoissenet, Waves called solttons (Springer, Berlin- 
Heidelberg, 1994). 



10 



Niurka R. Quintero, Angel Sanchez: ac driven sine-Gordon solitons: dynamics and stability 



2. A. Sanchez and A. R. Bishop, SIAM Review, in press 
(1998). 

3. Yu. S. Kivshar and B. A. Malomed, Rev. Mod. Phys. 61, 
763 (1989). 

4. A. Sanchez and L. Vazquez, Int. J. Mod. Phys. B 5, 2825 
(1991). 

5. Most important in contexts where soliton stability is a cru- 
cial factor, as it happens in many applications (Josephson 
junctions, fiber optics). Dynamical solitons are easy much 
more sensitive to perturbations; an archetypal example are 
those of the Korteweg-de Vries equation (see, e.g, |3| and 
references therein). 

6. G. L. Lamb, Rev. Mod. Phys. 43, 99 (1971). 

7. T. H. R. Skyrme, Proc. Roy. Soc. London A247, 260 
(1958). 

8. T. H. R. Skyrme, Proc. Roy. Soc. London A262, 237 
(1961). 

9. U. Enz, Phys. Rev. 131, 1392 (1963). 

10. R. Rajaraman, Solitons and Instantons (North-Holland 
Publishing Company, Amsterdam, 1982). 

11. A. Barone and G. Paterno, Physics and applications of the 
Josephson effect (Wiley, New York, 1982). 

12. E. Feldtkeller, Phys. Stat. Sol. 27, 161 (1968). 

13. J. D. Weeks and G. H. Gilmer, Adv. Chem. Phys. 40, 157 
(1979). 

14. A. Sanchez, D. Cai, N. Gr0nbech- Jensen, A. R. Bishop, 
and Z. J. Wang, Phys. Rev. B 51, 14 664 (1995). 

15. J. Frenkel and T. Kontorova, J. Phys. USSR, 1, 137 
(1939). 

16. F. R. N. Nabarro, Theory of crystal dislocations (Dover, 
New York, 1987). 

17. C. Cattuto and F. Marchesoni, Phys. Rev. Lett. 79, 5070 
(1997). 

18. S. W. Englander, N. R. Kallenbach, A. J. Heeger, J. A. 
Krumhansl and S. Litwin, Proc. Natn. Acad. Sci. USA 77, 
7222 (1980). 

19. M. Salerno, Phys. Rev. A 44, 5292 (1991). 

20. L. V. Yakushevich, Reviews of Biophysics 26, 201 (1993). 

21. L. V. Yakushevich, Studia Biophysica 121, 201 (1987). 

22. P. J. Pascual and L. Vazquez, Phys. Rev. B 32, 8305 
(1985). 

23. A. Sanchez, A. R. Bishop and F. Dommguez-Adame, Phys. 
Rev. E 49, 4603 (1994). 

24. D. Cai, A. R. Bishop and A. Sanchez, Phys. Rev. E 48, 2, 
1447 (1993). 

25. M. B. Fogel, S. E. Trullinger, A. R. Bishop and J. A. 
Krumhansl, Phys. Rev. B 15, 1578 (1977). 

26. D. W. McLaughlin and A. C. Scott, Phys. Rev. A 18, 1652 
(1978). 

27. P. S. Lomdahl, O. H. Olsen and M. R. Samuelsen, Phys. 
Rev. A 29, 350 (1984). 

28. J. C. Ariyasu and A. R. Bishop, Phys. Rev. B 35, 3207 
(1987). 

29. D. W. McLaughlin and E. A. Overman, Phys. Rev. A 26, 
3497 (1982). 

30. P. S. Lomdahl and M. R. Samuelsen, Phys. Rev. A. 34, 
664 (1986). 

31. P. S. Lomdahl and M. R. Samuelsen, Phys. Lett. A. 128, 
427 (1988). 

32. O. H. Olsen and M. R. Samuelsen, Phys. Rev. B 28, 210 
(1983). 

33. L. L. Bonilla and B. A. Malomed, Phys. Rev. B 43, 11 539 
(1991). 



34. D. Cai, A. Sanchez, A. R. Bishop, F. Falo and L. M. Fiona, 
Phys. Rev. B 50, 9652 (1992). 

35. P. J. Martinez, F. Falo, J. J. Mazo, L. M. Fiona and A. 
Sanchez, Phys. Rev. B 56, 87 (1997). 

36. N. R. Quintero and A. Sanchez, preprint cond-mat/- 
9709063, Phys. Lett. A, to appear. 

37. A. R. Bishop, D. W. McLaughlin and M. Salerno, Phys. 
Rev. A 40, 6463 (1989). 

38. W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. 
Flannery, Numerical Recipes in Fortran 2nd Edition (Cam- 
bridge University Press, Cambridge, 1992). 

39. For J?o = 0.9, the last point in Fig. pi the breather does de- 
cay into radiation instead of a into kink-antikink pair. We 
have included this point in the plot because in this case we 
could make sure that the breather had indeed disappeared. 

40. N. R. Quintero, A. Sanchez, and F. Dommguez-Adame, 
Resonances, bifurcations, and kink dc motion in the ac 
driven (f> 4 equation, to be published. 

41. J. C. Fernandez, M. J. Goupil, O. Legrand, and G. 
Reinisch, Phys. Rev. B 34, 6207 (1986). 



